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We calculate the general solutions for a warm inflationary scenario with weak dissipation, review- 
ing the dissipative dynamics of the two-fluid system, and calculate the bispectrum of the gravi- 
tational field fluctuations generated in the case where dissipation of the vacuum potential during 
inflation is the mechanism for structure formation, but is the sub-dominant effect in the dynamics 
of the scalar field during slow-roll. The bispectrum is non-zero because of the self-interaction of the 
scalar field. We compare the predictions with both those of standard, or 'supercooled', inflationary 
models, and warm inflation models with strong dissipation and consider the detectability of these 
levels of non-Gaussianity in the bispectrum of the cosmic microwave background. We find that the 
\£) , levels of non-Gaussianity for warm and supercooled inflation are an order of magnitude different. 
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^ I I. INTRODUCTION 

t-H ' 

The standard model of Hot Big Bang, or Friedmann-Robertson- Walker, cosmology with a simple, single-field 
supercooled inflationary phase, is highly supported by large scale structure 0], cosmic microwave background Q and 
weak lensing [j| observations, where supercooled inflation describes inflation without interactions between fields, ft 
is the most popular and developed model for use in parameter estimation, using the general assumption that the 
^ resulting perturbation spectrum is close-to scale free and Gaussian. The scenario has long been known, however, to 
\Q , contain a small but non-zero non-Gaussian effect due to the self-interaction of the inflaton field. This signal was 
0^ ' quantified in its contribution to the bispectrum statistic in Q, neglecting the effect of parametric amplification at 
inflation exit 0, . 

The challenge remains to make inflation, and cosmological parameter estimation using inflation, self-consistent, 
and to constrain parameters within every physically likely inflationary scenario. Warm Inflation Q describes the set 
• of models where the presence of interactions between fields during slow-roll is included in the analysis of inflation. 
Oh 1 The d\ mamics and observational consequences of warm inflation have been explored in several recent papers 0, 0> 
ITol ITlI Il2| . The production of radiation from the dissipation of the scalar field potential, which results from the 
presence of interactions, can lead to scenarios with levels of radiation density which are approximately stable. It has 
been shown that inflation can occur in the presence of a thermal component to the energy in the Universe. Several 
viable mechanisms for implementing such dissipation during inflation have been proposed in [T3I ITiL ITU IHI . The 
popular standard, or supercooled, inflation scenario is one of the limiting cases of warm inflation. In that case there 
is no radiation present either due to a lack of radiation production or due to the levels of radiation produced being 
redshifted away in the exponential expansion of the Universe during inflation. 

A second limiting case for the warm inflation scenario is inflation with strong dissipation where the interaction 
terms are greater in magnitude than the Hubble friction term of standard slow-roll inflation. In 0, El fTsl. fli^ and [jj 
solutions for the equation of motion of the scalar field have been calculated for the case of strong dissipation du ring 
inflation. There is a further limiting possibility out of this set of models, the weak dissipation scenario |2(1 1 12 lL |22( 
where a constant radiation component is created from dissipation which is the sub-dominant effect in the dynamics 
of the inflaton during slow-roll. 

In this paper we will present an analysis of the fluid dynamics of warm inflation which includes a more complete 
derivation of the energy constraints on warm inflation than has been previously outlined (i.e. in |20||). The solutions 
for the equations of motion of the scalar field and its spectrum of perturbations are then calculated for warm inflation 
with weak dissipation for polynomial potentials. In Section IIIII of the paper we use the solutions to provide an 
estimate of non-Gaussianity of the perturbations from warm inflation with weak dissipation which can be compared 
with the predictions for strong dissipation which are quantified in |9( for various warm inflationary potentials, and with 
the predictions for supercooled inflation calculated in Q . Consistency tests of the form of the fluctuations generated 
together with accurate predictions of non-Gaussianity should be an important tool in discriminating between structure 
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formation theories. We will conclude in Section with a discussion on the validity of these analyses in the wake of 
a recent preheating calculation for supercooled inflation j2^| which suggests that the non-Gaussianity from standard 
supercooled inflation has been previously greatly underestimated and may be of a level already ruled out by high 
resolution astrophysical data. 



II. THE BACKGROUND DYNAMICS OF THE 2-FLUID SYSTEM DURING WARM INFLATION 



In warm inflation, as is the case in supercooled inflation, the presence of the scalar field dominates the dynamics of 
the system. The energy density and energy pressure of the scalar field have the form: 
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Vffl (1) 



P4> = y - Vtf), (2) 

with V </> 2 /2. The potential energy of the scalar field dominates the system. The equation of motion of the scalar 
field can also be represented by the scalar field energy being divided into two fluids. The first fluid, representing the 
dynamical part of the scalar field, has an energy pressure equal to its energy density. The second fluid, representing 
the greater part of the energy, has an energy pressure equal to — lx its energy density. Interactions and dissipation 
of the scalar field during slow-roll are, in this analysis, modelled by a constant cross-section dissipation term. So the 
equation of motion of the scalar field during warm inflation has the form 

4> + 3H<j) + r<j) + V'{(j)) =0, (3) 

where the overdots represent time derivatives, H = is the Hubble parameter; R(t) is the cosmic expansion factor, 
r is a friction term which incorporates the interactions of the scalar field with an environment of particles, as warm 
inflation models inflation occurring in a heat bath (see e.g. |24j\ T does not represent the inflaton decay rate at the 
end of inflation in the supercooled case, which would be modelled by rapid oscillations of the inflaton field at the base 
of the inflaton potential. 

To analyse the evolution of the energies in the system as a result of the constant radiation production during 
this scenario of inflation we require a 2-fluid model. We label the two fluids in the following calculations with the 
subscripts a — {<f>, r} corresponding to the scalar field energy density and the radiation energy density. The stress- 
energy conservation equations of the scalar field and of the radiation must contain an energy transfer term, Q a , 
between the fluids. 

p a = -3H(p a +p a ) + Q a , (4) 
such that J2 a Q a = °- From Eq.0 and Eq.(0) 

P4> = -3H4> 2 + Q h (5) 



The equation of motion Eq.© substituted into Eq.JSJ together with the slow-roll assumption of inflation (j> <C 4> 
gives the form of the energy transfer term. 
Thus 

p r = -4Hp r + T4?. (7) 

The evolution equation of the inflaton reduces to 

alj, _ 1 dV(4>) 

dt m + r dd> u 
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FIG. 1: Evolution with efold of inflation, N, of scalar field f2^> and radiation fl r energy densities for Y / H = fO 



and the slow-roll condition to (3H + T)\<p\ ^> \(f>\. 

The dissipation causes radiation to be produced continuously from conversion of scalar field vacuum energy. In 
a preliminary work leading to the development of the warm inflation scenario poj . comparing temperature and 
thermalization scales with the Hubble scale during inflation it was shown that a dissipative component as small as 
r > 10~ 5 H was adequate to realise T > H. 

We have evolved the fluid interactions of Eas. (|5l7|l with a numerical simulation, where the dynamic part of the 
scalar field energy makes up 1% of the total inflaton energy. The figure demonstrates that once the system is in the 
state p r ~ 0, it remains in this state until the slow-roll conditions are violated. 

In this paper, expressions will be obtained which apply to the weak dissipative regime, r -C H , so that the dominant 
friction term in the equation of motion is due to the Hubble expansion, and Eq.JSJ) is further reduced to 



d(j) _ V 



(9) 



while the radiation temperature resulting from weak dissipation, T > H , means that we are dealing with a thermalized 
system in every efold of inflation. 

We can now evolve the behaviour of the inflaton in the weak dissipative warm inflation regime for a range of 
polynomial potentials. 

Focussing first on the background dynamics of the inflaton field in this regime, we start with this general form for 
the potential: 



(10) 



in the region < <f> < M where M is Mqut = 10 14 GeV. 
Solving Eq.© for these potentials leads to, for n ^ 2 



<j>{t) = M 



M"" 2 (n-2) Xt 
(n- 1)! 3~ff 



and for n 



(11) 



(f>(t) = M exp 



3H 



(12) 



III. THE EVOLUTION OF INFLATON PERTURBATIONS DURING WARM INFLATION 



We assume that the production of radiation during inflation will influence the seeds of density fluctuations, and we 
hold that this will apply when the temperature during inflation is greater than the Hubble parameter, T > H , and 
that in these bounds the thermal fluctuations of the scalar field must dominate over the quantum fluctuations. 
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In order to treat the fluctuations of the inflaton field 5(j)(x,t), we utilise the same stochastic approach which has 
been used in the calculations of non-Gaussianity which have been carried out for strong dissipative warm inflation 
and supercooled inflation 0, 0] . And also corresponding to the approach of and Q the effect of the self-interaction 
of the scalar field is treated preferentially above the back-reaction on the Hubble constant. It is assumed that the 
fluctuations are small — the full inflaton field is expressed as </>(x, t) = </>o (t) + <5</>(x, t) with <fio being the homogeneous 
coarse-grained averaged field, S(j)(x,t) <C (j>o(t), the scale over which the average is taken being larger than any scale 
of cosmological interest. From this the equation of motion for the full inflaton field with fluctuations emerges as 



d ^dT- = JH [ e ^ v2 ^M) - V'Mx,t)) + r,(x,t)] , (13) 



where r/(x, t) is a Gaussian noise term modelling thermal fluctuations with a mag nitude {k F T/2ir 2 )2 [ig| from the 
inflaton effectively sitting in a heat bath, kp is the wavenumber of a perturbation as it exits the inflationary causal 
horizon. 

The fluctuation-dissipation theorem determines the properties of the noise, as in With respect to physical 
coordinates and in momentum space, these properties, modelled in the simplest case are 

(n) = o (14) 



(f](k,t)f](k',t')) = 2rT(2 7 r) 3 ( 5 (3) (fc - k')S(t - t'). (15) 

The choice of physical coordinates arises since we are interested in the evolution of the inflaton mode while they 
are sub-horizon scale and the heat bath controls the fluctuations. 

In order to calculate the predicted bispectrum from Ea. (|13|) . the inflaton fluctuations must be expanded to first 
and second order: <5</>(x, t) = S(j>i(x,t) + <502(x, t), where 34>i — 0(S(j)) and (502 = 0(8<f) 2 ). 

We will follow the modes through their fluctuations inside the horizon, and through freeze-out at each scale of 
fluctuation. On this scale comparative energetics control the scalar perturbation. We treat the physical modes 
adiabatically with respect to the characteristic macroscopic time scale, the Hubble time ~ 1/H. Therefore it is 
sensible to use physical and not comoving momenta: 

^phys — ^comf- ■ (1^) 

The evolution of each mode can be integrated through steps of Hubble time ~ 1/H, and the combination of Hubble 
time-steps will result in the complete solution. Hereafter for physical momenta the notation will have no subscripts 
fcphys = k, and magnitudes will be denoted without boldfacing as k = \k\. 
Within each efold of inflation, inside the horizon: 



4(^1 («!,*)) = ^[-k 2 5^{k,t) - V"(^ (t))^i(fc,t) 
dt SH 

+ V (k,t)} (17) 

±-{5<p 2 {k, t)) = -L [-k 2 6Mk, t) - V"(Mt))5Mk, t) 
dt SH 

~lv"'(Mt))SMk,t) 2 }. (18) 

Dividing cosmic time into successive time intervals of order 1/H, t n — t n -\ — 1/H, the solutions of Eqs. (|17|) and 
(|TH|l for < n _i < t < t n are respectively 



<tyi(M) 



A(k, t - t n _ x ) [ dt /1 %^A(k, a - u^y 1 

+A(k, t - t n _ 1 )5Mke- H{t "- t "- l) ,t n _ 1 ) (19) 
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i4(fc,t-t„_i) / di 7 
dp 3 



B(t') 



(2tt) 3 



i(p,0^i(fe-p.O 



A(k,t'-t n -iY 



+A(k,t- tn-O^tfce-^*"-*-- 1 ),^-!) 



where 



j4(fc, t) = exp 



(20) 



(21) 



B(t) = 



V'"(Mt)) 
3H 



(22) 



IV. THE WARM INFLATION BISPECTRUM 



Single field inflation models broadly predict Gaussian primordial density fluctuations. Multiple field inflation models 
may lead to a non-Gaussian distribution (e.g. chi-squared) . When second-order effects are taken into account, 
however, there are corrections to these general predictions. There is a resulting non-Gaussian signal in the cosmic 
microwave back gro und radiation (CMB), the magnitude of which varies depending upon the self-interaction of the 
inflaton field HM El- 

A perfect Gaussian distribution would imply that the perturbations in the energy density have no connected 
correlations higher than the 2-point correlation function in real and Fourier space. 

We will now calculate the bispectrum statistic, the three-point correlation function of the density perturbation 
distribution in Fourier space. This quantity simply translates to the harmonic bispectrum of the CMB |27j . 

The leading order contribution to this three-point correlation function comes from two first-order and one second- 
order fluctuation as 



(6(f>{k 1 ,t)Scf>(k2,t)6(t>{k 3 ,t)) 



A(k 3 ,t- t 60 
dp 3 



1/H) 



Le>o 



A- l [kz,t' -t 60 -l/H)B(t') 



%0— 1/H 

{6fa(kuh)6<h.(p, O)Wi(fe,* 2 )<y0i(fc3 - p,t')} 



+A(k 3 ,t-t 60 - l/H)(S4>i(k 1> t 60 )6(l) 1 ( k 2,t6o) 

^(fcae-Steo-l/ff)) 

+(fei «-> fe 3 )) + (fea <-» k 3 )). 



Thus Eq. (|23J becomes 



(6cj>(ki,t)6ct>(k2,t)Scj>(k 3 ,t)} w B(t m )At F 
dp 3 



(2tt) 3 
(fei ^ fe 3 )) + [k 2 



i (fei , ti)5t/n (p, t')){8(j> 1 [k 2 ,t 2 )Hi (k 3 - P, 0) 



(23) 



(24) 



It is possible to write a general expression for the bispectrum for slow roll, single field, supercooled inflation models 
as well as for the set of warm inflation models in terms of products of power spectra and a prefactor, Ai n f, which 
depends upon the interactions during inflation, on the shape of the inflationary potential and upon the shape of the 
triangles used to measure the bispectrum statistic. 



(#(fci)$(fca)$(fc 8 )> = 

A inf (27r) 3 (S 3 (fci + fe 2 + fe 3 ) [P*(fei)P*(fe 2 ) + perms] 



(25) 
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The relation between the scalar field fluctuation and the gravitational field has the simple form (from |28|) 

*(fc)=-|^(fc), (26) 

5 (j) 

thus j4inf calculated using equilateral bispectrum triangles for a weakly dissipative warm inflation regime is 

v"i(Mt F )y 



a weak 




(27) 



which, plugging in the temperature T calculated using the relation p r = (g*7r 2 /30) T 4 , where g* is the number of 
relativistic fields ~ 150, and taking into account the relation between the radiation energy density and the scalar field 
potential ~ V(4>), given in Ea. (|l()(l . and the CMB amplitude, given by 

5 H = ~S<f>, (28) 

o (p 

which, from CMB results of the Wilkinson Microwave Anisotropy Probe, is 5h = 1-94 x 1CP 5 [2^, l30t l3l| . gives for 
a quartic potential a value of Aj^ ak = 7.65 x 1CP 3 . Our derivation of Ai n f mirrors the derivation of this value, the 
coefficient of the bispectrum of the scalar field gravitational potential, and also the coefficient of the bispectrum of 
the CMB in the Sachs-Wolfe region, of 0. 

V. CONCLUSIONS 

We have followed the behaviour of the inflaton field and of radiation during warm inflation in the limit of weak 
dissipation. We have evolved the scalar field perturbation up to second order and calculated the non-Gaussian 
signal in the bispectrum statistic which results from the self-interaction of the scalar field. The value we obtain 
^inf ak = 7. 65 x 10~ 3 is an order of magnitude smaller than the corresponding signal for a quartic inflationary potential 
for the strong dissipative and cool inflationary scenarios 0, |j| . 

At the time of the publication of 9] it was believed by the authors that the difference in the prediction for 
the bispectrum statistics of warm inflation with strong dissipation: A?i on8 — 7.44 x 10~ 2 ; and the prediction for 
supercooled inflation: A™ f o1 = 5.56 x 1CP 2 4]; would be overwhelmed by a contribution from general relativistic 
second order perturbation theory of Ai n { ~ 1 |32|. However, an energetic analysis of the evolution of second order 
perturbation theory [3j| has since shown that this is not the case. 

However, one incentive for the use of the warm inflation model over the supercooled inflation model, the implications 
of the preheating mechanism on the radiation perturbations, has a new possible observational implication. In | 23f it 
is demonstrated that parametric resonance during preheating can enhance the non-Gaussianity such that ^4?°f o1 can 
easily become as large as ~ 100. The implication of this for the result of this paper is that the difference between 
the predictions of Gaussianity for strongly and for weakly dissipative warm inflation may be both significant and 
observable in the future. Although this will not be distinguishable with the simple observation of the bispectrum with 
satellite CMB data, the best predicted constraint on this from the upcoming Planck satellite being Av4j n f = 10 [34| . 
the difference between the non-Gaussian signals from the strongly and weakly dissipative warm inflation scenarios 
might be tested for with others of the many non-Gaussianity statistics i.e. the trispectrum [35j and the two-point 
correlation of peaks ■ 
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